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ABSTRACT 


Far-infrared (FIR)-radio correlation is a well-established empirical connection 
between continuum radio and dust emission of star-forming galaxies, often used 
as a too l in determining star-formation rates. Here we expand the point made by 
Mur phy (120131) that in the case of some interacting star-forming galaxies there is 
a non-thermal emission from the gas bridge in between them, which might cause a 
dispersion in this correlation. Galactic interactions and mergers have been known 
to give rise to tidal shocks and disrupt morphologies especially in the smaller of 
the interacting components. Here we point out that these shocks can also heat 
the gas and dust and will inevitably accelerate particles and result in a tidal 
cosmic-ray population in addition to standard galactic cosmic rays in the galaxy 
itself. This would result in a non-thermal emission not only from the gas bridges 
of interacting systems, but from interacting galaxies as a whole in general. Thus 
both tidal heating and additional non-thermal radiation will obviously affect the 
FIR-radio correlation of these systems, the only question is how much. In this 
scenario the FIR-radio correlation is not stable in interacting galaxies, but rather 
evolves as the interaction/merger progresses. To test this hypothesis and probe 
the possible impact of tidal cosmic ray population we have analyzed a sample 
of 43 infrared bright star-forming interacting galaxies at different merger stages. 
We have found that their FIR-radio correlation parameter and radio emission 
spectral index vary noticeably over different merger stages and behave as it would 
be expected from our tidal-shock scenario. Important implications of departure 
of interacting galaxies from the FIR-radio correlation are discussed. 


Subject headings: c 
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INTRODUCTION 


It has long been known that there is a tight correlation between far-infrared (FIR) and 


(van der Kruit 

1971; 

Helou et al. 

1985; 

Condon 


19921; mm et ahl l200ll) . This correlation has been s hown to hold ov er almost five orders of 


magnitude for galaxies, not just at local redshifts flYun et al. 120011 ). but also for redshifts 


Ibar et al. 

2008; 

Ivison et al. 

2010; 

Sargent et al. 

2010) 


luminosities (e.g. 

physical reason behind the FIR-radio correlation is not fully und erstood (see e.g. IVoelklll989 


Helou fe Bicavl Il993t lLacki et al.l 120101; ISchleicher fc Beckl 120131) . It is thought to be due to 
the ongoing star-formation - dust absorbs light from the massive young stars and emits it in 
the FIR band, while galactic cosmic-ray (GCR) electrons accelerated in supernova remnants 
emit in radio band. Consequently, the FIR -radio correlat ion has been proven to be a powerful 
tool for determining star-formation rates (Condon 1992). 


D espite numerous studies that have claim ed that the FIR-radio correlation is relatively 
stable (ISargent et al.ll201(J: iBourne et al.ll20111) . there are also several contemporary observa¬ 
tions both at low and high-redshifts that question this. For example, some observations have 
shown that the tight linear FIR-radio correlation varie s in the case of galaxies in rich clusters 


( Murphy et ah 2009jRandriamampandrv et ah 2015 ). but also for distant starburst galax¬ 


ies 


(jSaiina et al. 20081) and distant sub-mm galaxies (S MGs), which were found to be radio 


brighter with respect to t he local FIR-radio correlation (IKovacs et al. 


2010l; Smolcic et al. 2014). Though it is generally considered (ISargent et al 


2006 


Magnelli et ah 


2010) that this 


correlation does not evolve with redshift, the value measured in the few samples where 
radio-loud active galactic nuclei have be en excluded was found to be more than 10% lower 
(IMagnelli et al.ll2010t ISaiina et al.1120081) . However, what is obvious is the scatter of data 
around this correlation. This scatter is thought to originate in part from 1) yo ung active 
galactic nuclei in which the radio activity has begun only recently (iDrake et al.l 12003 ). 2) 
from much stronger magnetic fields in starbursts than is suggested by the minimum energy 
estimate ([Thompson et al.ll2006l ). or 3) from ongoing interactions between galaxies (IMurphv 
2013h . 


The departure from the typical FIR-radio correlation and excess of radio emission 
was also found in the case of so-called ’’taff y” syste ms (in terac ting galaxies with strong 
synch rotron-emitting gas bridge between them; Condon et alil2002 ). As suggested by Murphy 
(1201 31) this excess in non-thermal emission is probably du e to part icle acceleration in large- 
scale shocks in bridges between interacting galaxies (Lis enf e ld &Volk 2010). I 11 this work 
we expand on this and point out that such departure from the typical FIR-radio correlation 
may not only be the case for bridges and ’’taffy” systems, but from cosmic rays accelerated 
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in tidal shocks in the galaxies themselves. 


Shock waves that arise during galactic interactions (mergers and close fly-bys) have 
been known to impact the star-formation history of interacting galaxies by t rigger ing star- 
formation and even leading to a starburst phase (see e.g. Sanders fe Mirabe l 1996 ). It was 
pointed out that tidal shocks that accompany these interactions can give rise to a population 
of tidal cosmic rays (TCRs). This can h ave a pote ntially significant impact on nucleosyn¬ 
thesis of light elements such as lithium (IProdanovic et al.l 2013|). Moreover, it is clear that 
the p resen ce of such cosmic-ray population will also result in an enhanced radio emission 
(Prodanovic et ah 2013). As a result, the FIR-radio correlation will be impacted. This may 
possible cause the dispersion seen in the relationship and have important implications for 
st ar- format ion measurement s. 


In this work we explore the effects that the presence of tidal shocks might have on the 
FIR-radio correlation in interacting star-forming galaxies and test this hypothesis on a small 
sample of interacting systems. 


2. FIR-RADIO CORRELATION IN INTERACTING GALAXIES 


The FIR-radio correlation is described using the qui parameter (Helpu et ah 1985) as 

Tftr \ , f S lA 


<?IR = log 


3.75 x 10 12 Wm- 2 


log 


Wm- 2 Hz- 1 


( 1 ) 


where S\ A is the continuum radio emission flux at 1.4 GHz per frequency such that S\ A oc v~ a 
and a is the radio spectral index, positive in vast majority of sources. Ffir is the rest-frame 
FIR dust emission flux. Yun et ah (j200lh analyzed the sample of 1800 IRAS galaxies and 
measured this value to be = 2.34 ± 0.01 with a dispersion of 0.25 dex. Both normal 
star-forming spirals and merger-induced luminous infrared galaxies were included in their 
sample. 


In this work we revisit the idea of galactic interactions being (in part) the source of 
the large dispersion of the relationship. As galactic interactions can give rise to tidal shocks 
in the interstellar medium of interacting galaxies, they can consequently produce a cosmic- 
ray population, in addition to normal GCRs, resulting in excess radio emission. Such an 
enhancement in radio flux would cause the dispersion of the FIR-radio correlation, which 
would especially be significant at higher redshifts with increased rate of galactic interactions 
compared to low-redshift Universe. 


Unlike GCRs, cosmic rays accelerated in tidal shocks due to close galactic fly-bys would 



















4 


result in an excess of non-thermal emission in gamma-rays and radio band that is not imme¬ 


diate ly accompanied by the corresponding increase in star-formation rate (Prodanovic et al. 


2013). Therefore the most promising way for identifying the presence of TCRs is to ob¬ 


serve galaxies in early (mid) merging stages (merging stage 3 according to Haan et al. 2 011 


classification scheme). One would also anticipate additional heating of the gas and espe- 
cially dust in tidal shocks (a s was observed in interacting systems; see e.g. lLanz et al.ll2013l; 


Mentuch C oo per et ah 20.12), compared to what would normally be expected due to the 


ongoing star-formation, again, provided that the system is observed in an sufficiently early 
merging stage. With these additional sources of non-thermal and thermal radiation it is clear 
that the FIR-radio correlation would be affected such that giR ^ q]- R - a typical FIR-radio 
parameter will not be equal to that same parameter qJ R for a system where there is a presence 
of additional cosmic-ray population (such are for example TCRs). More specifically, if there 
are tidal shocks and a TCR population present, then the expected parameter would be 


T 

9ir 


log 


FIR 


rTS 

WlR 


C I cTS 
‘-’1.4 ~r *5 1 4 


( 2 ) 


The Fj Ri is the additional FIR flux of dust coming from tidal shock heating, while Sjf 
is the additional radio flux from TCR electrons. Observing that giR > 1 and assuming 
that the effects of tidal shocking are a small perturbation already existing effects (i.e. that 
^fui/SJa < F Fm /S 1A ) we will have that in general g IR > qf R . However, as the interaction 
between galaxies progresses, how would we expect this parameter to change? 


2.1. Early Interaction: Enhanced Heating 


At the very early stages of interaction, we expect tidal shocks to form in the ISM and 
start heating the dust and gas. Dust can be heated in collisions with gas or by shock 
UV radiation, causing it to radiate thermally in infrared. However, this will also result in 
destruction of dust due to sputtering processes. The time-scale of destruction of dust in 


shocks goes from few thousand years to tens of milli ons of years (ITielensI Il998t lLau et al. 


20151 : IVillar-Martm et al.l l200ll: iDwek fe Arendtl 119921 1 . depending mainly on the grain size, 
strength of the shock and density of the ambient medium. On the other hand typical 
cosmic-ray acceleration timescale in supernova remnants is of the order of lifetime of the 


remnant, being a few ~ 10 5 yr (e.g. Strong et all 120071 ). For large-scale shocks considered 


here, cosmic-rays will be accelerated as long as the tidal shock propagates, which can be of 
the order of Gyrs. Thus, during the first few thousand to million years of galactic interaction, 
we can consider that there is enhanced thermal emission of heated dust and gas over what 
is typically expected of non-interacting systems, without enhanced non-thermal emission of 
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freshly accelerated cosmic-rays. As a result, we can expect a non-thermal radio emission 
dominated by the already present GCR population, and thermal emission with an added 
contribution from tidal shock heating, leading to q ir < g IF [ . The enhanced overall spectral 
index will also ch ange, becoming shallower (harder) a < 0.8, compared to a typical observed 
value a ~ 0.8 (jCondonl 119921). A shallower radio sp ectral index indicates a more dominate 
thermal component (see e.g. L isenfeld feV olk 2000i ). and vice versa. We expect this change 
in spectral index to accompany the change in the FIR-radio parameter on a same timescales. 

It is generally accepted that galaxies like the Milky Way are not to be considered as 
closed-box systems (unlike starburst gala xies), but should be expected to ’’leak” cosmic- 
rays on timescales of r esc ~ 2 x 10 7 yr flGarcia-Munoz et al.lll977lh So for a g alaxy like 
the Milky Way, with supernova rate of Rsn ~ 1/50 yr (Tannnann et al. 1994), 4 x 10 5 
supernova events will occur before escape losses become important, allowing an estimate of 
the maximal GCR flux. In order for tidal shocks to result in a cosmic-ray flux greater than 
the already present GCR flux, their input needs to be equivalent to about 10 5 supernova 
events. Consequently, the volume of gas shocked by tidal shocks needs to be equal to volume 
shocked by that many supernova events, assuming the same acceleration efficiency. Taking 
tha t parti cle s are efficiently accelerated in supernova remnants up to the radius of Rsn ~ 10 
pc (Berezhko & Volk 1 2004 1. we find that tidal shocks would need to shock all the gas up to 
the radius of about 1 kpc in a galactic disk of thickness d = 300pc. This would take about 
10 million years, assuming a conservative radial tidal shocks with velocity u s hock ~ lOOkm/s 
(tidal shocks are weaker than collision sh ocks w hich are of theonier of the velocity of merging 
galaxies ~ few hundred km/s; see e.g. Kashi yarn a & Meszaros 20.14). Therefore, we estimate 
that this early stage where non-thermal emission due to TCRs can be neglected, would last 
during the first few million of years, or more precisely 

3/2 / , \ - 1/2 

loji 1 ^shock \ / 1 i-KN \ / & 

Tearly £ 13Myr 


f -Rsn V 


X 


T, 


hshock 

100km/s J lylOpc ) 

1/2 / p 

/ -fts 


300pc 


i-SN 


2 x 10 7 yr J ll/50yr 


(3) 


In the case of starburst galaxies with ~ 10 x higher supernova rates, this timescale 
would be about a factor of 3 higher. Note however that supernova (i.e. star-formation) 
rate determined from radio observations would not be suitable within the framework of our 
hypothesis because if tidal cosmic-ray population is present it might lead to enhanced radio 
emission resulting in a overestimate of star-formation rates. A better way, free from possible 
” contamination”, would be to use star-formation rates determined from Ha observations to 
test how galactic interactions might impact our understanding of star-formations rates. This 
will be the topic of a follow-up study. 
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2.2. Mid Interaction: Enhanced Particle Acceleration 


At few tens of million of years since the beginning of interaction, acceleration of particles 
in tidal shocks will start to be significant. In an extreme scenario, non-thermal emission of a 
population of tidal cosmic-rays will at some point become dominant over the galactic cosmic- 
ray population. This would be late enough in the interaction so that the star-formation has 
been triggered by tidal shocks, resulting in additional heating of the gas and dust due to new 
stars. The approximate start of this stage would be at the timescale of stellar contraction 
and formation, that is, this stage would begin at least tens of millions of years after the 
formation of tidal shocks. 


Tidal shock heating then becomes less important as the source of thermal emission. 
On the other hand, TCRs are now a dominant source of non-thermal emission given that 
formation of new stars has not yet been accompanied with the increase in GCRs flux. In that 
case we expect a decrease in the FIR-radio parameter below its typical value gj^ cr < giR, 
such that gj^ h > giR > qJ R > g/j/ cr . The decrease in the FIR-radio parameter would be 
accompanied with a steeper (softer) radio spectral index a > 0.8, compared to a typical 
observed value in radio, reflecting a larger presence of a non-thermal component. Again, this 
evolution in spectral index would be expected to be on the same timescale as the evolution 
of the FIR-radio parameter. 


A quick look into energetics also leads to the conclusion that tidal cosmic-ray population 
may become dominant at some point. Namely, one supernova event injects approximately 
EgN = 10 51 /50 = 2 x 10 49 erg/yr into the interstellar medium, where ~ 10% of that energy is 
generally considered to go into particle acceleration. We note that this energy-injection rate 
is valid for normal, star-forming galaxy like the Milky Way. In the case of starburst galaxies 
with much higher supernova rates, energy-injection rate would clearly also be higher. This 
means that the following estimate is valid for the early epoch of galactic interaction, where 
supernovae from subsequently triggered starburst phase have not yet injected their energy. 
That will be the epoch where tidal cosmic-ray flux can be a pote ntially sign ifica nt addition to 
the already present galactic cosmic-ray flux. Following Prodano vic et ah 0201311 we estimate 
the kinetic energy of the encounter between Milky Way type galaxies at a distance of 50 
kpc is E mSL j = 4 x 10 6 ° erg. For a minor merger where smaller component is 1000 times less 
massive energy of interaction is E min = 4 x 10 57 erg. Thi s energy estimate is of the order 
of tidal energy th at get s injected in the system (see e.g. Spitzer 1958). If 10% (see e.g. 


Kash i vama & Meszaro s 2014) of this energy goes into particle acceleration and if interaction 


timescale is of the order of 10 9 yr, we find that the rate at which tidal shocks can inject 
energy into particles is fhr ima j = 4 x 10 5 ° erg/yr for major merger and fhr im in = 4 x 10 4 ' 
erg/yr for minor merger. We see that in the case of major mergers, tidal shocks have the 
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potential to dominate the non-thermal emission. In the case of minor mergers TCRs can be 
at the 10% level of the GCRs, which is a smaller effect. Note however that this was estimated 
assuming that the smaller component has large star-formation and supernova rates equal to 
the Milky Way. In reality, for a smaller system, we would expect this rate to be orders of 
magnitude lower for a galaxy 1000 times less massive than Milky Way. In that case energy 
injection rate of tidal shocks compared to supernovae would be much higher. On the other 
hand, in later phases, when star-formation is triggered by the merger, star-formation rate of 
such small systems can indeed be as high as for the Milky Way. 


Another way to look into energetics is to compare the energy of tidal interaction with 
the non-thermal luminosity of the source. Let us, for example, consider the case of stage 3 
system NGC5256, with component masse s of 6.3 x 1O 1O M 0 and 5.3 x 1O 1O M 0 and separation 
of about 7.5 kpc (IMazzarella et ah 2012). If we then assume that they tidally interact over 
the timescale of 1 Gyr, we find that tidal energy injection rate is of the order of 7% = 
3 x 10 42 er g/s. If we c ompare this to source luminosity at 4.8 GHz vL v = 4.1 x 10 39 erg/s 
(IMazzarella et ahl 2012) we see that for electron-to-proton ratio of 1/100 at energy 1 GeV, 
~ 10% of tidal interaction energy needs to be converted into 1 GeV cosmic-rays to account 
for luminosity at 4.8 GHz band. Acceleration efficiency for even this cons ervat i ve est im ate 
is consistent with typical acceleration efficiencies in supernova remnants (Lisenfeld &Volk 

201flh . 


In order to check if TCRs accelerated in large-scale tidal shocks that traverse through 
entire interacting galaxies, can produce observed non-thermal radio emission and spectra, 
we analyze the case of IC1623, which is one of the systems in our sample that is a stage 3 
merger with lower than typical FIR-radio parameter qui = 2.08, and steeper than typical 
spect ral index a = 0.91. Following Bell’s theory of particle acceleration in shock fronts (iBe ll 
19781 ) we calculate radio emissivity of shocked gas assuming that the entire gas of IC1623 was 
shocked by tidal shocks of speed 250-300 km/s, that number den sity of I SM is 1 cm -3 and 


that the average galactic magnetic held is of strength 10-25 /j,G (jDrzazga et ah 20111; Beck 


2009 ). Based on that, and adopting its measured spectral index we calculate the expected 
luminosity at 1.4 GHz L 14 = 1.3 — 4 x 10 23 W/Hz. This is in very good agr eement with 
luminosity calculated from observed flux L li4 = 2.1 — 4 x 10 23 W/Hz (iMurphy 120131 ). 


2.3. Late Interaction: Enhanced Star-Formation 

Tens to hundreds of million of years from the beginning of interaction, star-formation 
rates would be increasing, accompanied by the increase in supernova rate but with slight 
offset at order of lifetime of massive stars. This is consistent with results from numerical 





















models (e.g. Di Matteo et al. 2008) which give that first enhancements in star-formation 
rates would be few hundred million years in the interaction after the first passage, followed 
by the star-burst epoch at time ~ 1 Gyr from the beginning of interaction and lasting some 
500 Myrs. During that time the FIR-radio parameter would start to grow as FIR emission 
would be increasing due to new stars being born. Finally, at the end of interaction one would 
expect a similar epoch of enhanced supernova rate. At that point, the FIR-radio parameter 
would be expected to go back to its typical value comparable to that of isolated systems. 
Finally, we note that at higher redshifts, timescales and durations of each p hase would be 


different due to lower concentration of dust and lower metallicity (e.g. Pettini et al. 1997) 


Table 1: Merger stage classification according to Haan et ah (2011) with corresponding sym¬ 
bols used on plots. 


Merger 

stage 

Symbol 

Description 

0 

None 

Non-merger 

1 

A 

Pre-merger: separate 
galaxies, non tidal tails 

2 

A 

Ongoing merger, early: 
progenitor galaxies distinguishable 

3 

• 

Ongoing merger: 
progenitors sharing an envelope 

4 

* 

Ongoing merger, late: one galaxy 
with double nuclei and a tidal tail 

5 

□ 

Post-merger: one galaxy with single 
(disturbed) nucleus and prominent tails 

6 

None 

Post-merger, late: one galaxy 
with single nucleus and weak tail 


3. Analysis of The Sample of Interacting Galaxies 


As described, what would eventually be expected from galactic interactions is the de- 
pendance of the FIR-radio parameter g IR on the merger stage. Early in the interaction this 
parameter would have larger than the average value for the non-interacting star-forming 
galaxies, after which it would be expected to decrease to lower than the average value, fol¬ 
lowing eventual relaxation back to the original value. To test t his we have analy zed the 
sample of 43 IR bright galaxies found in data sets from Dopita et ah (2002) and Murphy 
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(2013). D ata presente d in 


Murphy (2013) were drawn from IRAS revised Bright Galaxy 


Sample (Sanders et ah 2003). Galaxies in this sample were chosen to have 60/nn flux densi- 
ties larger then 5.24 Jy and FIR luminosities > lO 11 ' 25 !/^. Taking into account that origina l 


Dopita et al.l (120021) data set consists of two studies (IKewlev et al.ll2001t ICorbett et al.ll2002l) 
where different lower limit for FIR luminosities was chosen, our sample used in this work 
satisfies the overall criteria that all chosen galaxies are IR bright with 60/um flux densities 
larger than 2.5 Jy and have IR luminosities > 1O 1O ' 5 L 0 . We have excluded all sources clas¬ 
sified as active galactic nuclei in the original data sets. All galaxies presented here have well 
sampled radio spectra between 1.4 and 8.4 GHz, but for most of them dat a is also available 
outside of that range of frequencies. In total 15 ga laxies from iMurphvI ( 20131) are used. 
Values of giR were taken directly from Murphy (120131 ). For Dopita et al. (2002) sample we 
h ave calculated Qm following equation (1), estimating the far-infrared flux, Fpim as defined 


m 


Sanders & Mirabel (jl996i) . based on IRAS flux densities f v at 60/im and 100/irn: 


Wm 


-2 


FIR ' = 1.26 x 10~ 14 


2.58/60um + /lOO 


J/xm 


Jy 


(4) 


For galaxies from Murphy ([2013) radio spectral indices were recalculated for archival 

4.8/1.4). Galaxies 


VLA data at 1.4 and 4 . 8 GH z from standard a = log (Sj a / Sa a) / log 


studied in IDopita et ahl (2002) are COLA galaxies (Corbett et ah 2002), and their indices 


were determined at same frequenci es as in the original 


Corbett et al. 


stages have bee n defined follow ing lHaan et ahl (120llh. Classification status for each galaxy 


was taken from Murphy 02013 ) and 


(120021) paper. Merger 


Dopit a et al. (120021 ). We have excluded galaxies unclas- 


QP 1 

sided by their merging stage in iMurphvl ( 20131 ). while for one of them (NGC 6286S) we have 
additionally confirmed (inspecting the optical DSS images on the NASA/IPAC Extragalac- 
tic Database - NED) its merging stage 3 according to apparent weak tidal features and an 
envelope shared with its interacting companion. 

We have also calculated the mean values and standard deviations of gm parameter and 
spectral index a separately for all interaction stages of galaxies in our sample. This is 
presented in Table 2. On Figure [j] we plot FIR-radio parameter g IR as a function of the 
merger stage for this sample of IR bright interacting galaxies. Th e origina l i dea that the 
FIR-radio parameter qm is sensitive to merger stage was presented in Mu rphy (12013! ). Figure 
5 from Mur phy (120131 ) revealed significant scatter of the FIR-radio parameter for sources at 
different merger stages, and here, we have improved the statistics with additional data set 
of IR bright sources to see is the same behavior of qm related to merger events will hoick 
Though statistics is still limited, we see that only in the case of stage 4 and 5 mergers is 
the_mean value of the FIR-radio parameter consistent with its global value of giR = 2.34 
(Yun et al. 2001). The average giR for the other merger stages deviate from the global value. 
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Merger stage 


Fig. 1.— Plotted is the FIR-radio correlation parameter gjR for the sample of 43 interact¬ 
ing star-forming galaxy pairs at different merger stages. Merger stages fjHaan e_t_ ah 2011) 
were labeled using the sam e sym bols as noted in Table 1. The dashed line represents the 
typical value of the giR flYun et all 2 0011 ). while dotted lines represent its mean deviation. 
Mean values and their standard deviations were calculated for each merger stage subsample 
separately and are presented with red symbols and blue error bars. 
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For example, for the most significant case of merger stage 3 (where we expect that the 
impact of tidal cosmic-ray population is most prominent), we perform a statistical t-test 
to check how significant the offset of this data subset is is from the sample mean value of 
< 7 ir = 2.34 ± 0.21. We find a t-value of t — 3.8 corresponding to a p-value of p — 0.004, 
indicating that the probability that this offset is random is less than 1%. This offset is even 
more significant if we compare it to the population mean with gj R = 2.34 ± 0.01. What it 
appears is that systems in close fly-bys, in the pre-merger stages, have higher than average 
values of q \ R , which then decreases toward later merger stages and reach es the minimum 
value at merger 3 systems, after which it increases again. We note that Murphy (2013) has 
shown similar trend plotting the difference between the observed and nominal logarithmic 
IR and radio flux densities, but without any further statistical analysis. We should point 
out that two of merger stage 3 galaxies show ’'no excess” to standard gi R value. However, 
inspecting NED images of those sources reveals that source IRAS 06295 shows no obvious 
tidal features, which is a typical property of stage 3 mergers. The other galaxy IRAS 12286- 
2600 has a smaller companion with double nuclei, which is again not typical of a simple stage 
3 system. 


As was discussed in Section 2, the evolution of the FIR-radio parameter g IR with respect 
to the merger stage would also be accompanied with the corresponding evolution in the radio 
spectral index. Ongoing mid-stage mergers would be expected to have enhanced non-thermal 
emission with ax > a = 0.8 and < q = 2.3. This trend is actually seen on Figure [2] where 
we have plotted radio spectral index a as a function of the FIR-radio parameter gi R . The 
upper left quadrant of the plot where gi R is lower than typical a = 0.8 ± 0.05 (dashed line) 
and radio spectral index is higher than typical (dashed line), is dominated by merger stage 3 
and 4 systems. As opposed to that, the bottom right quadrant with gi R higher than typical 
and radio spectral index lower than typical value, is dominated by the late merger stages 
4 and 5, as would be expected. In the case of merger stage 3 the t-test gives a p-value of 
p = 0.03 that this offset of stage 3 spectral index from typical value is random. 


Table 2: FIR-radio parameter gi R and spectral index a values calculated for each merger-stage 
subsample separately. 


Merger 

Symbol 

<?IR 

a 

1 

A 

2.47 ± 0.12 

0.69 ± 0.171 

2 

O 

2.38 ± 0.03 

0.56 ± 0.185 

3 

• 

2.09 ± 0.21 

0.92 ± 0.19 

4 

* 

2.31 ± 0.17 

0.79 ± 0.22 

5 

□ 

2.44 ± 0.15 

0.51 ± 0.16 
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1 
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Fig. 2.— Plotted is the radio spectral index a (Si. 4 oc u~ a ) measured between 1.4 GHz and 
4.8 GHz, as a function of the parameter for the same sample of interacting star-forming 
galaxies at different merger stages (Haan et ah 20111 labeled using the same symbols as 
noted in Table 1. Das hed lin es represe nt typical values of radio spectral index a = 0.8 
flCondon 1992) and qm (|Yun et ah 20 01 ) parameter for star-forming galaxies. For our data 
set we have calculate the scatter around the_mean spectral index to be a = 0.8 ± 0.05 which 
is consistent with results from Murphy (120130 . while the scatter of the mean FIR-radio 
parameter is gm = 2.34 ± 0.21. 
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4. DISCUSSION AND CONCLUSION 


In this paper we explore and draw attention to important effects that close-galactic 
interaction and mergers can have on the stability of the widely-used FIR-radio correlation. 
Close galactic interactions produce tidal shocks in interacting galaxies which leads to gas 
and dust heating. Moreover, particle acceleration in tidal shocks gives rise to a tidal cosmic- 
ray population. These two effects would impact the infrared and radio emission of the 
interacting galaxies and cause variations in the FIR-radio parameter q\R and radio spectral 
index over different merger stages. To test this, we have analyzed the sample of 43 IR bright 
interacting galaxies in different merger stages, looking at how their FIR-radio parameter 
and radio spectral index change over different merger stages. What we have tentatively 
found is that Qir first decreases during early merger stages, and then later increases. This 
is consistent with what would be expected if the heating of the gas and dust in tidal shocks 
is the dominant effect at early merger stages, followed by the phase where tidal cosmic- 
ray emission dominates over the existing GCR emission, and eventually ending with the 
enhanced star-formation taking this correlation parameter back to its typical value. As a 
result, if interacting galaxies are included in the FIR-radio correlation, they could be one 
possible cause of its observed dispersion, and probably a dominant cause of its dispersion 
going to a higher redshifts where merger rates are higher. 


Our results are consistent with results of Mur phy (120131 who was the first to p oint 
out the possible impact of galactic interactions on the FIR-radio relation. While Mu rphy 
(120131 ) was mostly focused on describing the radio emission of ” Taffy”-like systems, here we 
wanted to also include the cases where there is amplified radio emission like in the galaxy 
NGC5256 and galaxies alike which cannot be treated as ” Taffy”-like due to lack of H 2 gas 
in the bridge between the two galaxies. More spe cifically, the NGC5256 consists of a pair 
of galaxies, comparable in size and scale (IPetrosian et al.1119801 ) with projected separation 
between their nuclei of ab o ut 8 kpc. It has a very low qm = 1.90 value, and recent multi- 
wavelength studies (IMazzarella et ah, 120121 uncovered several interesting features in these 
luminous IR-galaxies: (a) The optical morphology of the NGC 5256 north-east nuclear 
environment is similar to the radiative shock observed south of the nucleus of M51; (b) 
Bridge of CO(l — 0) emission is spatially decoupled from the radio continuum emission; (c) 
Steeper radio index is observed not only in the region in betwee n galaxies, but also close 


to the edges of both components (see radio maps in Vardoulaki et al. 2 Ola ): (d) Whereas 


the bulk of the HCO + molecules in Taffy systems is located between the nuclei, the HCO + 
in NGC5256 is still bound to the galaxies; (e) Spectral energy distribution (SED) of NGC 
5256 can be modeled as a starburst-dominated, however compared to galaxies with the same 
SED description (for example, NGC 2623 or NGC 6240), but different merging stage, the 
NGC5256 shows evidence of higher dust temperatures; (f) a soft X-ray emission extending 
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15_kpc to the north of the system, between the nuclei, was observed (Brassington et ah 


20071 ). revealing the presence of shock-heated gas indicating that corresponding synchrotron 
radio emission is shock-induced. This case shows that even in non-” Taffy’” systems there 
can be important sho cking of gas and dust and particle acceleration within the interacting 
galaxy. Furthermore, Drzazga et al. (2 0110 have studied how tidal interactions affect the 
evolution of galactic magnetic fields and according to their analysis of polarization data 
in two Taffy systems (UGC12914/5 and UGC813/6), they found a well-ordered magnetic 
held in the bridge in one of them (UGC12914/5), corresponding more to a s h ocked gas we ll 
before a adiabatic (Sedov) phase, which was the scenario analyzed in Lisenfeld feVolkl (12010! ) 
to explain the bridge emission. All of this indicates that tidal shocks in interacting galaxies, 
their evolution and effects can provide a general solution to effects observed in interacting 
and merging systems. 


If tidal shocks and cosmic-ray acceleration are the underlying reason for merging and 
interacting galaxies to deviate from the well-established FIR-radio correlation, this effect 
can be used at high redshifts as a tool in searching for interacting systems and testing our 
understanding of high redshift interaction rates. Moreover, departure from the FIR-radio 
correlation in the case of inter acting systems could have important consequences for determi 


natio n of star-formation rates (mm fc Carillijl2002c lBellll2003l: ICaxilli et a.l .1120081: IPunne et ah 


2009 ) leading to its overestimate. A more reliable way to determine star-formation rates 


m 


interacting galaxies would be to look into Ha emission that directly probes emission of 
massive young stars. For example, the mean star-formation rate determined from radio ob¬ 
servations for our entire sample of 43 galaxies is 41M 0 /yr, while Ha observations give the 
average of 16M 0 /yr. The effects of this will be explored in the follow-up work. 


In order to determine if indeed interacting systems deviate from the well-established 
FIR-radio correlation and how, a larger sample of bright IR galaxies along with radio data 
has to be analyzed. The data available from surveys such are COSMOS and CANDELS 
along with the upcoming data from ALMA will provide a perfect testing ground for this. 
For example, CANDELS survey is using deep near-IR imaging to reveal morphological clas¬ 
sifications of distant IR-bright galaxies and directly counts the number of interacting pairs 
up to (z ~ 2), while COSMOS survey uses interferometric follow-up observations of distant 
dusty star-forming galaxies to redshifts even greater than z ~ 4. With ALMA it would 
be possible to resolve these dusty galaxies into individual pairs for a large sample of radio 
faint submillimeter galaxies (SMGs). Combining different radio maps it would be possible to 
determine spectral indices of distant SMGs, and trace the FIR-radio correlation for a large 
sample of SMGs free from the biases. 


Besides the need for a larger sample, it would also be important to obtain multi- 
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wavelength observations of systems that are found to be good candidates to be dominated 
by tidal cosmic ray population. An example of such systems would be the Whirlpool galaxy, 
M51 (o r N GC 5194) and NGC 5256. Recent PdBl Arcsecond Whirlpool Survey (PAWS) 
(Schinnerer et al. 2013|) revealed the presence of additional cosmic ray emission in spiral 
arms of the NGC 5194 where low rate of ongoing star formation is present. The smaller 
companion, NGC 5195, was found to have low or no ongoing star formation flBigiel et ah 


2008 ), however, it has uncharacteristically high dust temperature (jMentuch Cooper et al. 


2012]), possibly due to shock heating. 


Since the effects of tidal shocks would be most pronounced in the smaller of the interact¬ 
ing components for which very little data is currently available, it would also be important 
to observe best candidates of those smaller galaxies especially in radio and 1R domains. 
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